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The Loop-Cluster Algorithm for the Case of the 6 Vertex Model* 
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We present the loop algorithm, a new type of cluster algorithm that we recently introduced for the F model. 
Using the framework of Kandel and Domany, we show how to generalize the algorithm to the arrow flip symmetric 
6 vertex model. We propose the principle of least possible freezing as the guide to choosing the values of free 
parameters in the algorithm. Finally, we briefly discuss the application of our algorithm to simulations of quantum 
spin systems. In particular, all necessary information is provided for the simulation of spin i Heisenberg and xxz 
models. 
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1. INTRODUCTION 

Cluster algorithms, originally introduced for 
the Ising model 0| and then generalized to vari- 
ous other situations |^,|| , are one of the promising 
ways of overcoming critical slowing down. Re- 
cently PI PI we introduced algorithms for vertex 
models MM , which are the first cluster algorithms 
for models with constraints. While Q is an adap- 
tation of an algorithm originally devised for solid- 
on-solid models, the loop algorithm introduced in 
js) does not resemble any existing scheme. 

In H] we presented the loop algorithm for the 
F model. This enabled us to present our idea as 
clearly as possible. Here we shall show how to 
generalize it for the 6 vertex model, which has an 
additional coupling and a richer phase structure 
(see below). The framework of |^] proves to be 
an extremely useful tool here. 

Our scheme is devised such as to take into ac- 
count the constraints automatically, and to al- 
low a simple way to construct the clusters. After 
defining the relevant probabilities, we find that we 
still have some free parameters. In order to op- 
timize the algorithm, we introduce the principle 
of minimal freezing. As seen from the example 
of the F model, this choice of parameters is of 
utmost importance. 

The loop algorithm can be further generalized 
to more complicated vertex models. However, as 



it stands, there are already important applica- 
tions: quantum spin systems can be simulated by 
mapping them to vertex models In particu- 
lar, the loop algorithm for the 6 vertex model pre- 
sented here can be used to simulate spin ^ Heisen- 
berg ferromagnets and antiferromagnets, and xxz 
models, even in more than one dimension. 

2. THE 6 VERTEX MODEL 



The six vertex model [|6|,[7| is defined on a square 
lattice. On the bonds there lives an Ising-like 
variable that is usually represented as an arrow. 
For example, arrow up or right means plus one, 
arrow down or left means minus one. At each 
vertex, there are two incoming and two outgoing 
arrows. In fig. 1 we show the six possible config- 
urations at a vertex, numbered as in [Hj7|. 

The statistical weight of a configuration is given 
by the product over all vertices of the vertex 
weights p(u). For each vertex there are 6 pos- 
sible weights p{u), u — 1,...,6. We assume the 
vertex weights to be symmetric under reversal of 
all arrows. So in standard notation MM we have: 



p(l) = p(2) = a , 
p(3) = p(4) = b , 
p{5) = p(6) = c . 



(1) 



*To appear in Lattice '92, Amsterdam 1992, ed. J. 
et al., Nucl. Phys. B (Proc. Suppl.). 
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The six vertex model basically has two types 
of phase transitions: of Kosterlitz-Thouless type 
and of KDP type A sub- model exhibiting 

the former is the F model, defined by c = 1, a ~ 



1 2 3 4 5 6 

Figure 1. The six vertex configurations, u = 1, ...,6 (using the standard conventions of |^). 



h = exp(— i<r), K > Q. For the latter transition 
an example is the KDP model itself, defined by 
a = 1, 6 = c = exp (-if), K>Q. 

3. THE LOOP ALGORITHM 

If we regard the arrows on bonds as a vec- 
tor field, the constraint at the vertices is a zero- 
divergence condition. Therefore every config- 
uration change can be obtained as a sequence 
of loop-flips. By "loop" we denote an ori- 
ented, closed, non-branching (but possibly self- 
intersecting) path of bonds, such that all arrows 
along the path point in the direction of the path. 
A loop-flip reverses the direction of all arrows 
along the loop. 

Our cluster algorithm performs precisely such 
operations, with appropriate probabilities. It 
constructs closed paths consisting of one or sev- 
eral loops without common bonds. All loops in 
this path are flipped together. 

We shall construct the path iteratively, follow- 
ing the direction of the arrows. Let the bond b 
be the latest addition to the path. The arrow on 
h points to a new vertex v. There are two out- 
going arrows at v, and what we need is a unique 
prescription for continuing the path through v. 
This is provided by a break-up of the vertex v. 
In addition to the break-up, we have to allow for 
freezing of v. By choosing suitable probabilities 
for break-up and freezing we shall satisfy detailed 
balance. 

The break-up operation is defined by splitting v 
into two pieces, as shown in fig. 2. The two pieces 
are either two corners or two straight lines. On 
each piece, one of the arrows points towards v, 
while the other one points away from v. Thus 
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Figure 2. The three break-ups of a vertex: U-ur 
(lower-left-upper-right), ul-lr (upper-left-lower- 
right), and straight. 



we will not allow e.g. the ul-lr break-up for a 
vertex in the configuration 3. If we break up v, 
the possible new configurations are obtained by 
flipping (i.e. reversing both arrows of) the two 
pieces independently. On the other hand, if we 
freeze v, the only possible configuration change is 
to flip all four arrows. 

The break-up and freeze probabilities are con- 
veniently described within the general framework 
for cluster algorithms proposed by Kandel and 
Domany It is sufficient to give them for one 
vertex, which is in the current configuration u. 
We define 6 new interactions (weight functions) 
Pi, i = 1, 6, corresponding to specific break-up 
and freeze operations. (The labelling of the new 
interactions is completely arbitrary, and the fact 
that we have six of them is just a coincidence). 
For each vertex in configuration u, we replace 
with probability Pi{u) the original interaction p 
by the new interaction pi. Detailed balance and 
the proper normalization of probabilities require 
that for every u 



p{u) ^ 



(2) 
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action 


1 

freeze 1,2 


2 

freeze 3,4 


3 

freeze 5,6 


4 

U-ur 


5 

ul-lr 


6 

straight 


Pt{u) 


1, p{u)=a 
0, else 


1, p{u) = h 
0, else 


1, p{u)^c 
0, else 


0, p{u)^a 

1, else 


0, piu) = b 

1, else 


0, p{u) = c 

1, else 




qi/a, p{u)^a 
0, else 


92/6, p{u)^b 
0, else 


93/c, p(m) = c 
0, else 


0, p{u)^a 
qi/p{u), else 


0, p{u) = b 
qs/piu), else 


0, p{u)^c 
q(i/p{u), else 



Table 1. New interactions Pi{u) and the probabilities Pi{u) to choose them at a vertex in current con- 
figuration u. See eq. (||). 



where > are parameters. 

As discussed in |^] (see also table 1), freezing 
is described by introducing one new interaction 
for each different value of p{u). For example, to 
freeze the value a, we choose the interaction pi 
to be pi{u) = 1 if p{u) = a, and pi{u) — other- 
wise. In other words, if m is 1 or 2, the Boltzmann 
weight pi{u) is one, so transitions between 1 and 
2 cost nothing; the vertex configurations 3, 4, 5, 
and 6 are however not allowed with pi. 

Each break-up is also described by one new in- 
teraction. As an example take the ul-lr break-up. 
It is given by the new interaction number five, 
with P5{u) = 1 if p{u) = a or c, and P5{u) — 
if p{u) = b. In other words, with the new inter- 
action p5, transitions between 1, 2, 5 and 6 cost 
nothing, while the vertex configurations 3 and 4 
are not allowed. This corresponds precisely to al- 
lowing independent corner fiips in a ul-lr break- 
up (see figs. 1,2). 

The full list of new weights Pi{u) and proba- 
bilities Pi(u) to choose them are given in table 1. 
From (H) we also obtain: 

91 + 95 + 96 = o , 

92 + 94 + 96 = , (3) 

93 + 94 + 95 = c . 

Assume now that we have broken or frozen all 
vertices. Starting from a bond 69, we proceed 
to construct a closed path by moving in the ar- 
row direction. As we move from vertex to ver- 
tex, we always have a unique way to continue the 
path. At broken vertices the path enters the ver- 
tex through one bond and leaves it through an- 
other. If the last bond b added to the cluster 
points to a frozen vertex v, the path bifurcates 
in the directions of the two outgoing arrows of 



V. One of these directions can be considered as 
belonging to the loop we came from, the other 
one as belonging to a new loop. Since we also 
have to fiip the second incoming arrow of v, we 
are assured that this new loop also closes. The 
two loops have to be flipped together. In general, 
the zero-divergence condition guarantees that all 
loops will eventually close. 

We have now finished describing the procedure 
for constructing clusters. In order to specify the 
algorithm completely, we must choose values for 
the constants qi, and decide how the clusters are 
flipped. The former problem is of utmost impor- 
tance, and it is the object of the next chapter. For 
the cluster flips, we may use both the Swendsen- 
Wang procedure and the single cluster flip ||^ . In 
1^ we used the latter, and obtained a drastic re- 
duction of critical slowing down. 

4. OPTIMIZATION 

We have seen that freezing forces loops to be 
flipped together. Previous experience with clus- 
ter algorithms |^ suggests that it is advantageous 
to be able to flip them independently. We there- 
fore introduce the principle of minimal freezing 
as a guide for choosing the constants qi: we shall 
minimize the freezing probabilities, given the con- 
straints (H) and 9i > 0. In [|| we report that for 
the case of the F model, optimization by mini- 
mal freezing does indeed minimize critical slow- 
ing down. Here we discuss optimization for the 4 
phases of the 6 vertex model, usually denoted by 
capital roman numerals 

Let us first look at phase IV, where c > a + b. 
To minimize the freezing of weight c we have to 
minimize q^. From (^), q^ = c-a-b+qi+q2+2qe. 
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With Qi > this imphes (?3,min — c — a — b. 
The minimal value of qs can only be chosen if 
at the same time we set qi — q2 — 0, i.e. mini- 
mize (in this case do not allow for) the freezing 
of the smaller weights a and b. The optimized 
parameters for phase IV are then: 

91 = 0, 92 = 0, qs = c - a - b, 
(?4 = fe, 95 = a, 96 = . 

In phase I the situation is technically similar. 
Here a > 5 + c, and we minimize freezing with 
qi = a — b — c and 92 = 93 = 0. The same holds 
for phase 11, b > a + c, where we obtain minimal 
freezing lor q2 = b — a — c and 9i = 93 = 0. 

Phase III (the massless phase) is characterized 
by a, 6, c < ^{a + b + c). Here we can set all 
freezing probabilities to zero. Thus, 

91 = 0, 2q4 = b + c - a , 

92 = 0, 2^5 = c + a - 6 , (5) 

93 = 0, 2qe = a + b- c . 

The F model is obtained from (^ and (^ as 
the special case a — b < 1, c = 1. One can easily 
see that for this case we recover the discussion of 
Q. (Notice that since a = 6, in the F model the 
straight break-up will be called freezing). 

5. APPLICATIONS, CONCLUSIONS 

We have presented a new type of cluster algo- 
rithm, the loop algorithm, for the case of the six 
vertex model. For the F model, the algorithm has 
been shown in j^] to beat critical slowing down. 

Particularly promising is the possibility of ac- 
celerating Quantum Monte Carlo simulations 
1^. Quantum spin systems in one and two di- 
mensions can be mapped into vertex models in 
1 + 1 and 2 + 1 dimensions via the Trotter for- 
mula and suitable splittings of the Hamiltonian 
Q . The simplest example is the spin i xxz quan- 
tum chain, which is mapped directly into the 6- 
vertex model. For higher spins, more complicated 
vertex models result (e.g. 19-vertex model for spin 
one). 

For (2 + 1) dimensions, different splittings of 
the Hamiltonian can lead to geometrically quite 
different situations [|[||. We can e.g. choose be- 
tween 6-vertex models on a complicated 2 + 1 



dimensional lattice, and models on a bcc lattice 
with 8 bonds (and a large number of configura- 
tions) per vertex. Notice that for the simulation 
of the 2-dimensional Heisenberg antiferromagnet 
using the former splitting, all relevant formulas 
have been worked out in the present paper. 
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